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ABSTRACT 

We present X-ray observations of the afterglow of GRB000926 performed 
around and after the break observed in the optical light curve two days after 
the burst. The steep X-ray light curve observed around the break confirms 
the presence of this feature in X-rays. However, the spectral and temporal 
properties are not consistent with a standard jet scenario based on synchrotron 
emission, requiring a more complicated model. We find that X-ray and optical 
data are compatible with a moderately coUimated fireball (with opening angle 
9 ^ 25°) expanding in a dense medium {n ^ A x 10^ cm~^). This produces 
two breaks in the light curve. The first, at t ~ 2 days, is due to jet behavior. 
The second, around 5 days, is attributed to the transition of the fireball to a 
non-relativistic expansion. This transition predicts ^flattening of the light curve, 
that explains the late X-ray measurement in excess above the extrapolation of 
the simple jet scenario, and is also consistent with optical data. 

Subject headings: gamma rays: bursts 
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1. Introduction 

Afterglows of Gamma Ray Bursts (GRB) are typically characterized by a power law 
behavior in time and in frequency, i.e. the observed flux F ^ (t — tg)^^^ This property 
is nicely accounted for by the fireball model, in which the afterglow is produced by the 
interaction of a relativistic expanding shell with an external medium (e.g. IWijers, Rees fc| 



Meszaros 1997| ). The spectral and temporal slopes depend on the geometry (e.g. spherical 
vs jet ( [Sari, Piran fc Halpern 1999| )) and on the density distribution of the external medium 
(e.g. constant density - hereafter ISM - vs. wind ( [Chevalier fc Li 1999D ). A change of the 
slope (i.e. a break) of the light curve is produced either when a spectral break (like that 
associated to Uc, the cooling frequency of electrons ) transits the observed frequency window 
or when "bulk" variations,- i.e. global changes of the kinematics or hydrodynamics of the 
fireball,- take place. The latter is the case of jet expansion or the transition from relativistic 
to non- relativistic expansion (hereafter named NRE). In contrast to the effect of spectral 
breaks, jet expansion and NRE produce achromatic breaks in the light curves. In a jet this 
happens when the beaming angle of the fireball (~ F"^, where F is the Lorentz factor), 
becomes similar to the opening angle of the jet, 6 ([Sari, Piran fc Halpern 1999[ [Roads[ 
[T999[) , at the time te: 

to ^{l + zK^y^H^f 'days, (1) 
n 0.1 

where n = lO^cm"^ is the density and Ei = -Ej^ss 10^^ erg s~^ is the energy the 
fireball would have if it were spherically symmetric; i.e., in the case of a jet of solid angle 
An, its energy would he E = ^Ei. 

The transition to NRE instead occurs when the blast wave has swept up a rest-mass 
energy equal to its initial energy, at a time (e.g. jWijers, Rees fc Meszaros 1997[ , Livio & 
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Waxman (2000), hereafter LW): 



W^3i±^(^)^/='days (2) 
2 ne 

Both jet and NRE can therefore yield achromatic breaks a few days after the burst. 
The NRE requires an external medium of a much higher density than that of a jet. This is 
indeed what one would expect if GRB originate in dense star forming regions, a scenario for 
which supporting evidence is accumulating (e.g. [Piro et at. 20004 |Kulkarni et al. 2000| ). 



So far, evidence of achromatic breaks has been claimed in some optical afterglows 
(e.g. [Harrison et al. 1999| , p<.ulkarni et al. 1999| ), but an unambiguous detection of this 



feature in the X-ray range is still missing (e.g. [Kuulkers et al. 2000| , |Pian et al. 2001| ). 



Most of the authors attribute this feature to jet behavior, although this model may have 
some problem in accounting for the sharpness of the transition ([Rhoads fc l:'ruchter 2001] ). 
In fact, there is no - a priori - reason to exclude a NRE and there is at least one case 
(GRB990123) in which it has been attributed to either cases ( [Kulkarni et al. 1999 , Dai fc 



[Lu 1999|) . In GRB000926 a break in the optical light curve at t ~ 2 days has been reported 
by several groups ([Fynbo et al. 2000b| , |Veillet 2000| , [Halpern et al. 200q , |Rol et al. 2000| , 



Price et al. (2000, hereafter PO). The observed variation of the temporal slope of ~ 2 is 
much larger than the value of 0.25 expected by a spectral break ( [Sari, Piran fc Halpern 



1999| ). Indeed, a jet scenario has been proposed to explain this behavior (P0|). 



In this paper we report on BeppoSAX and Chandra observations of the afterglow of 
this burst (Sect. 2 & 3). Those observation occurred around and after the break observed in 
the optical. We compare X-ray with optical data (Sect. 4) and show that they are explained 
by jet behavior and a NRE into a high density medium. 
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2. Observations, data analysis and the X-ray spectrum 

The gamma-ray burst GRB000926 was detected by the Interplanetary Network (IPN) 
on to = Sep. 26. 99 U.T. ( [Hurley et al. 2000|) . It was a moderately bright GRB, with a 



revised 25-100 keV fluence of 6.2 x 10 ^erg cm ^, and a duration of about 25 s ( [Hurley 



\et al. 2000|) . Optical observations ([Pall et al. 2000] , [Gorosabel et al. 2000| ) showed an 



optical counterpart, at a very bright level {R = 19.3) one day after the GRB. This is the 
second brightest afterglow ever observed in the optical band. This finding triggered two 
rapid target of opportunity observations in X-rays, one with BeppoSAX and the other with 
Chandra. A second Chandra observation was scheduled around t = to + 12 days, to study 
the late time evolution of the afterglow. 



The BeppoSAX ( [Piro, Scarsi fc Butler l9U3 , Boella et al. 1997 ) observation started 



on Sep. 29.03 (t — to =2 days) and lasted for 12 hrs. Effective exposure times were 20 
ksec for the MECS(1.6-10 keV) and 5 ksec for the LECS(0.1-10 keV). The analysis of the 
MEGS image showed a previously unknown fading source in a position consistent, within 
the error of 50", with the optical transient, and therefore identified as the X-ray afterglow 
of GRB000926 ([Piro et al. 2000b| , [Piro et al. 2000^ ). The source was only marginally 



detected in the LEGS. The average count rate of the source in the MEGS(1.6-10 keV) is 
(3.2 ± 0.6) X 10~^cts/s. The spectrum is well fitted by a power law with energy index 
a = (0.9 ± 0.7)0 and absorption by our Galaxy (A^^^ = 2.7 x lO^^cm'^) and F(1.6-10 
keV)=F(0.2-5 keV)= (2.6 ±0.6) x lO^^^erg cm~^ s~^. We have checked different background 
subtraction methods, finding that the corresponding systematic effect amount to less than 
about 15%. This has been included in the flux error quoted above. 



^ Hereafter errors and upper limits on fit parameters, including temporal slopes, 
correspond to 90% confidence level for one parameter of interest 
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A Chandra X-ray Observatory ( |Weisskopf et al. 1996|) TOO observation (hereafter CI) 
with the ACIS-S began on Sep. 29.674 (t — 1^=2.1 days) and ended on Sep. 29.851, with 
an exposure time of lOks. The X-ray afterglow was observed at a position RA(2000)=17h 
04m 09.6s; Decl. (2000)= +51 47' 8.6" with an uncertainty of 2", consistent with the 
optical position. The source count rate in the 0.3-5 keV range is (2.9 ± 0.2) x 10^^ cts 
s^^, corresponding to a flux F(0.2-5 keV)= (1.2 ± 1.0) x lO^^^erg cm^^ s~^ for the best 
fit power law spectrum. This has a = (0.8 ± 0.2), i.e. consistent with the BeppoSAX 
slope, and absorption Nh = (5 ± 3) x 10^°cm~^, consistent with that of our Galaxy. 
Since the BeppoSAX and Chandra spectral shapes are consistent, we have fitted them 
simultaneously linking the spectral parameters (but for the normalizations). We obtain a 
good fit (x^ = 25.6, u = 26) with a = 0.85 ± 0.15, a ratio of the BeppoSAX over Chandra 
normalization of 2.3 ± 0.6 and some marginal evidence of intrinsic absorption (in addition 
to the one in our Galaxy) in the rest frame of the GRB (z=2.06 [Fynbo et al. 200^ , [Castro 



et al. 2000|) ), corresponding to a column density Nhgrb = {O-AILqII) x lO^^cm ^. The 



residuals around 2 keV, i.e. the iron line region in the rest frame of the GRB (Fig.l ), are 
marginal, with an upper limit on an iron line of / < 2 x 10~^cm~^ s^^, corresponding to 
an equivalent width EW < 1 keV. For comparison, the iron line detected by Chandra in 
GRB991216 has EW = 0.5 keV ( piro et al. 2"000i|) . 

The last Chandra observation (hereafter C2) took place on Oct. 10. 176 and ended on Oct. 
10.58. The afterglow count rate had faded to 2.0 ± 0.75 x lO^^cs"^ (0.3-5.0 keV), allowing 
~ 60 counts to be collected in the 33 ksec exposure. With so few counts it is difficult to place 
good constraints on Nh and a simultaneously, so we fix Nhgrb = (or 0.4) x lO^^cm"^ 
and find a = 0.85 ± 0.15 (or 1.2 ± 0.3)). The fiux in the 0.2-5 keV range derived by the best 
fit power law is equal in both cases to (8.3 ± 1.1) x lO^^^erg cm~^ s~^. An overall fit on the 
three data sets (see Fig. P gives x^/^ = 30.8/24, with spectral parameters consistent with 
those derived previously (see the inset in Fig.l) and a ratio of the normalizations C2/C1 of 
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(7.7 ±1.1) X 10-2. 

We find marginal evidence of a steeper spectrum in C2. Tlie ratio of count rates 
in the 0.2-1.5 keV range over 1.5-8 keV is in fact 2.1 ± 0.6 and 3.9 ± 0.2 in CI and C2 
respectively. If we leave the slope of C2 unlinked in the combined fit of all data sets, we get 
ac2 = 1-23 ± 0.3. While this result is of marginal statistical significance, it is interesting to 
note that the difference in slope between C2 and Cl-fBeppoSAX is Aa = 0.4 ± 0.3, i.e. 
consistent with the steepening of 1/2 expected when z/c passes into the energy window. 

An intrinsic column density Nhgrb ~ 0.4 x lO^^cm-^ is marginally detected at a 
confidence level of ~ 98% (see the inset in Fig.l). We remark that, since most of the local 
gas is ionized by GRB photons QBoettcher et al. 1999|) , the real column density can be 
much higher. Assuming that the relation between the optical extinction Ay and Nh is 
similar to that in our own galaxy (e.g. Predehl fc Schmitt 19951 ), gives ^ 2. This value 



is ~ 2 — 20 greater than that derived by optical measurements (P0|), implying a dust-to-gas 
ratio much lower than that of our Galaxy. This effect has been observed in other GRB 
( palama fc Wijers 2000| ) and it may be the result of destruction of dust by the GRB 



radiation ( |Waxman fc Draine 2000| ). 



3. A very steep X-ray decay 

In the BeppoSAX data the source exhibits a substantial decay, decreasing by a factor 
of (1.7 ± 0.5) in 6 hours. This corresponds to a power law decay F oc {t — to) ''^ with slope 
6x ~ 4. This behaviour is consistent with the first Chandra observation, that took place 
near the end of the BeppoSAX observation, and whose flux was a factor 2.3 ± 0.6 times 
lower than the average BeppoSAX flux. Combining these two data sets (Fig. 2) we derive a 
slope Sx = 3.7lJ;g, in agreement with the optical slope measured after the break. The late 
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Chandra point appears to deviate from this law, the slope connecting CI with C2 being 
5x = 1.70 ± 0.16. However, when the large error of the earlier slope is considered, the two 
measurements are nearly compatible (Fig. |]). A combined fit to all the X-ray data points 
gives then 5x = 1.89j:o;J6 with F(2-10 keV)=1.19 x lO^^^g^g cm'^ s"^ at t = 1 day. 

The X-ray decay of GRB000926 is much steeper than that typically observed in X-ray 
afterglows, where 6x ~ 1.4 QPiro 2000|) . We remark, however, that previous observations 



covered the evolution of the afterglow primarily at t — to < 2 days. The data presented here 
are instead sampling the X-ray behaviour at later times. It is in this domain that we expect 
deviations due to either a jet or NRE. 



4. Jet vs NRE of a spherical fireball in a dense medium 

The optical data (Fig.2) are nicely fitted by the standard jet model with synchrotron 
emission (|P0| ) but, when X-ray data are considered, we find inconsistencies with this 
scenario. The relationship linking the temporal and spectral slopes are (|Sari, Piran fc 
Halpern 1999|) 6 = 2a + 1 [y < Vc) and 5 = 2a = [y > Vc)- We find then that the X-ray 
temporal and spectral slopes 5x = 1.9 and ax = 0.85 are consistent only with the case of a 
cooling frequency below the X-ray range. The slope p of the electron distribution Nf, oc 
would then he p = 6 = 2a = 1.7 ± 0.3, i.e. marginally compatible with the minimum slope 
needed to avoid an infinite energy content in electrons. A more relevant problem of this 
model is that the predicted temporal slope is significantly lower than the one measured 
in optical (pO| ; see also Fig. |l]), 6o = 2.3 ± 0.08. Furthermore, the electron slope derived 
from the optical data under the same scenario {p = 2.39 ± 0.15) is not consistent with the 
X-ray derivation. Finally we note that the - admittedly marginal - evidence of a spectral 
steepening in C2, argues against a jet behavior, during which i>c remains constant. 
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We now examine the alternative explanation of a transition to NRE of a spherical 
fireball . The relationships linking 6 and a are given in Livio & Waxman (2000). We find 
that the X-ray data are consistent with a NRE either in a uniform density medium ( which 
gives 5 = 2.17 ± 0.45) or in a wind {6 = 2.3 ± 0.35). The cooling frequency is now above the 
X-ray range, and therefore p = 2.7 ± 0.3. In this case the optical and X-ray spectral data 
points should lie on a single power law. We have evaluated the slope apx = ~ ^iog{Jl//^o) 
the power law connecting the flux at 1 keV observed by Chandra {Fx) with the R band 
flux extrapolated at the time of the Chandra observation (Fo, Fig. 2) and corrected for the 
extinction (pO|). We find aox = 0.9 — 1.1, in agreement with the slope measured in X-rays. 
Before the break the evolution is indistinguishable from that of a spherical relativistic 
fireball, which follows a decay with 5 = 1.3 ± 0.2 (ISM) or 5 = 1.8 ± 0.2 (wind). The 
latter case is not compatible with the optical slope measured before the break (pO|), and we 
conclude that both optical and X-ray data are consistent with a spherical expansion in a 
uniform medium, with the temporal break produced by a transition to NRE. From eq.(2) 
we derive rig ~ 10. We have estimated, following Freedman and Waxman (2000), a fireball 
energy comparable to the observed 7-ray isotropic energy, which is Ei^^ = 8 x lO^^erg 
{z = 2.03 ( ICastro et al. 2000| ), Hq = 75 km/s/Mpc, go = 1/2). The absorbing column 
density is Nh ~ nrj^RE ~ 10^^ cm^^, where the radius tnre ~ lO^'^^cm ( pi;W|) . This is 
in apparent disagreement with the upper limits derived by the X-ray spectrum. This 
disagreement is rectified by recalling that the medium is heavily ionized by the GRB 
photons, and therefore its effective absorption is neglig ible QBoettcher et al. 1999|) . 



5. A collimated fireball in a dense medium 

There remains, however, a discrepancy between optical and X-ray data, that is not 
accounted for by previous models. The X-ray decay after the break {5x = 1-70 ± 0.16) is 
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significantly flatter than the optical one {So = 2.30 ± 0.08, FO). In fact, the extrapolation 



of the X-ray flux with the optical slope at the time of the second Chandra observation 
underestimates the observed flux by a factor of 2 ± 0.5 (Fig.|^). This indicates the presence 
of a flattening taking place in between the two sets of observations, as predicted in the case 
of a jet expansion followed by a transition to NRE ([LW|) . In this case, the lateral expansion 
of the jet begins when it is relativistic, and during the lateral expansion stage there is rapid 
decline of the flux due to the rapid decrease of the flreball energy per unit solid angle (as 
this angle increases). This is producing the flrst break in the light curve. The time at 
which the jet completes its side-ways expansion, is roughly the same time at which the flow 
becomes sub-relativistic. At this stage, lateral expansion is no longer important, and the 
rapid flux decline due to lateral expansion ceases, leading to a flattening of the light curve. 
This effect appears more evident in X-rays, which are not contaminated at later times by 
the presence of a fairly bright host galaxy (pO|). 

We show now that both optical and X-ray data are in agreement with this scenario. 
The light curve of the afterglow is composed by three subsequent power laws with decay 
slopes 6i, 62, and ^3, that we have modeled with the empirical formula : 



F{t) = F{tbi) ^ fjl\52s]-l/sll^(^\53s]l/s (3) 

The parameter s describes the sharpness of the breaks. The temporal slopes are related 
to the power law index of electrons by relationships that depend on the location of and 
on the density distribution of the external medium ( |LW|) . We can exclude most of the cases 
as follows, aox ~ «z implies that Uc is either below the optical or above the X-ray range. 
In the flrst case, we would have p = 2ax = 1-7 that, besides being less than 2, would 
give temporal slopes flatter than observed. In fact: 5i = ^^^Ez^ = q §2 = p = 1.7 and 
^3 = 1 + I(p - 2) = 0.65 (wind) or ^3 = 1 + |(p - 2) = 0.55 (ISM). During the jet phase 
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z/c = const while, in the NRE phase, i>c increases (decreases) with time for a wind (ISM). 
In the case of a wind Uc could move above the optical range, but still giving a flat slope 
'^3 = i + i(p-2) = 1.15. 

Let us now examine the case of lying above the X-ray range. This implies 
p = 2ax + 1 = 2.7. The case of a wind is excluded because 6i = 1.8, i.e. greater than 
observed. For a ISM we have 6i = |(p-l) = 1.3, 62 = p = 2.7 and ^3 = ^ + |(p-2) = 1.95, 
which appear to be consistent with the data. We have assessed this case more quantitatively 
by fitting eq.(3), with the slopes expressed as a function of p, to the optical data (with the 
contribution of the host galaxy subtracted) . We obtain p = 2.57 ± 0.3, in nice agreement 
with the value derived from the X-ray spectra, tbi = 1.76!to;4 days, tb2 = 5.5 ± 2 days, and 
s = 1.7 ± 0.3. Moreover, this same law (i.e. rescaling the normalization but with the same 
parameters determined by the optical light curve) fits fairly well all the X-ray data points, 
and in particular reproduce C2 (Fig. We also note that this law predicts a temporal 
slope of about 2.6 around 2 days, nearer to the steep decay observed in the BeppoSAX an 
CI observations. 

The simultaneous determination of = t^i and t^RE = h2 allows an unambiguous 
determination (i.e. not depending on estimates of Ei and n) of the jet angle 
( [TWD . In fact tNRE = (1 + zf'\^Y'K'^ which, when combined with eq.(l), 
gives 6 = 0.7( ^J^^ )^/^ = 0.4. The total 7-ray energy is then, for a two-sided jet, 
= O.lGE^i = 1.2 X lO^^g^g s-\ Finally, the density uq = (1 + zfEi^?,J^ ^ 0.04. 

' ''NRE 

Note that the jet angle is rather insensitive to the errors in the determination of break 
times. Conversely, the density should be considered as an order of magnitude estimate, 
due to the strong dependence on t^RE and to model details (|LW|) . We have finally checked 
whether the value of density is compatible with the broad band spectrum of the GRB, 
following the prescriptions of Wijers & Galama (1999). Since the parameters one derives 
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are dependent on the model adopted, we have restricted the analysis to the data before the 
break, when the fireball evolution is indistinguishable from the spherical case. Prom the 
radio data (Frail, private communication; Harrison et al. (2001)) and the extrapolation 
of the optical-to- X-ray spectrum at lower frequencies, we derive that the synchrotron 
self- absorption frequency Va and frequency of the maximum flux should be in the range 
between 10^° — lO^^Hz, and that the flux at z/^ should be ~ a few mJy. Furthermore, as 
we have shown before, 10^® Hz. The large range spanned by these quantities gives 
a density spread over orders of magnitude that, however, comfortably include the density 
estimate derived before. 

6. Conclusions 

In this paper we have shown that X-ray observations of the afterglow of GRB000926 
around and after the break observed in the optical confirm the presence of this feature in 
X-rays as well. However, they are not consistent with the standard jet scenario based on 
synchrotron emission, requiring a more complicated model. We find that both the spectral 
and temporal properties of X-ray and optical data are compatible with a moderately 
coUimated fireball {9 ~ 25°) expanding in a dense medium (n ~ 4 x 10^ cm~^). This 
produces two breaks in the light curve. The first is due, as usual, to a jet behavior. The 
second, around 5 days, is attributed to a transition to NRE due to the high density of the 
external medium, and is followed by a flattening of the curve, that accounts for the late 
X-ray measurement, and is also preferred by optical data. Optical measurements at later 
times should help to characterize this behavior but, to do so, it will be necessary to actually 
single out the afterglow from the host in the image. Current estimates of the contribution 
of the host galaxy from the flattening exhibited by the hght curve will, in fact, cancel out 
any intrinsic flattening due to the afterglow. Finally, we note that we cannot exclude other 
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possibilities like, for example, an event of gravitational lensing ( Mao fc Loeb 2001 ) that 



amplified the flux around 12 days, or a contribution of an Inverse Compton component 
to the late X-ray emission (Harrison et al. 2001). An analysis of different options goes, 
however, beyond the scope of this paper. 

Our result supports the association of GRB with star forming regions. In fact, the 
density is consistent with that of very dense molecular clouds {n ^ 10^ cm~^). On the 
other hand it is not as high as that required to produce the iron lines observed in some 
X-ray afterglows (e.g. |Antonelli et al. 200l| ; |Piro et al. 2000a] and references therein), which 



is possibly connected with much denser {n ^ 10^°cm~'^) progenitor ejecta. In these cases 
the environment must be composed by at least two regions, since the fireball requires a 
medium with a density n ^ lO^cm^^ to produce an afterglow lasting at least a few days. 
The properties of this medium, in particular the density, are not the same in all GRB. 
While in this GRB we find n ~ 4 x lO^cm"^, in several other GRB the density should be 
much lower, as indicated, for example, by the absence of breaks in the optical light curves 
( [Kulkarni et al. 2000| ) or by the evidence of relativistic expansion several weeks after the 
GRB ( [Waxman, Kulkarni fc Frail 1998 ). Another element of diversity in GRB appears to 



be the magnetic field strength. In this burst the magnetic field density is estimated to be 
less than ^ 10^^ times the equipartition value, comparable to that of GRB971214 (|Wijers ^ 



Galama 1999| ) and GRB990123 ( |Galama et al. 1999| ), but four orders of magnitude weaker 



than in GRB970508 QWijers fc Galama 1999| ). As suggested by Galama et al. (1999), such 
differences in field strength may reflect differences in energy flow from the central engine. 
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Fig. 1.— The X-ray spectrum of the afterglow of GRB000926 as observed by the ACIS-S 
of Chandra (first observation: open squares; second observation: filled squares ) and the 
MECS of BeppoSAX (filled circles; the LEGS is not shown for clarity). The continuous lines 
represent the best fit power law convoluted through instrument responses. Deviations from 
the best fit in a are plotted in the lower panel. In the inset we show the the contour plot 
of the photon index (a + 1) vs Nh with the 99% (thick hne) , 90% (medium-thick hue) and 
68% (thin line) confidence regions 
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Fig. 2.— The X-ray and optical (R band) light curves of the afterglow of GRB000926. 
Optical data are from Ball et al. (2000), Fynbo et al. (2000a, 2000b), Halpern et al. (2000), 
Price et al. (2000), Veillet et al. (2000). The X-ray light curves have been derived by 
converting the flux at 1 keV, which is not significantly dependent on the spectral shape, to 
the F(2-10) keV flux, by adopting the best fit power law with a — 0.85 for all the spectra. 
The dashed line is the best fit to the data of the BeppoSAX and first Chandra observations. 
The allowed region (at 90% confidence level) obtained by extrapolating the previous fit is 
identified with the shaded area. The continuous hue is the best fit to the optical data 
including the effect of the host galaxy (PO). The dash-dot hne is the template of the optical 
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Fig. 3. — Same as in Fig.2 but with the contribution of the galaxy subtracted from optical 
data. The continuous and dotted light curves represent the expected behaviour in the case 
of a moderately coUimated jet with a transition to non-relativistic expansion at t=5.5 days. 
Note the flattening of hght curves both in the optical and X-rays after this second break. 



